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' Abstract 

. We look at the possibility that appreciable neutrino masses and flavor mix- 

C*~) \ ing occur only within material media, driven by an interaction between leptons 

and a very light scalar particle. Limits are placed on the scalar particle mass 

. and coupling constants from a number of experimental and astrophysical con- 

! siderations. 
OO . 

^ . 1. Introduction. The recent results on atmospheric neutrinos [0-0] offer a per- 

Q\ \ suasive argument that neutrino mass and flavor oscillations are needed to resolve 

the anomalies with respect to standard model predictions. Once one accepts this 
conclusion, a parallel explanation of the solar neutrino deficits seems inevitable. Nev- 
ertheless, the evidence that neutrinos in a vacuum have mass or oscillate is still 
Qh! inconclusive. In particular, in the vacuum oscillation model that has been use to fit 

the data in 0, except for the predictions for the lowest energy bin, the flavor mixture 
of neutrinos from above is changed little from that at the neutrino production point, 
while the flavor mixture of upcoming //„ is changed greatly, by virtue of the much 
| greater path lengths from the point of production. In the lowest energy bin, p < .4 

GeV/c , the directional determination appears to be too imprecise to test the premise 
that only the upmoving neutrinos are changed. It is therefore worth considering the 
possibility that the effects depend not just on the distance but on the intervening 
matter as well. Indeed, Wolfenstein |J long ago discussed such a possibility, pointing 
out that experiments do not rule out flavor changing neutral currents, at essentially 
standard strength, that could cause neutrino flavor to precess in matter while not 
producing, or demanding, neutrino mass. 

In somewhat the same spirit we here suggest that the intervening matter induces 
effects that are tantamount to masses and flavor oscillations of the upmoving at- 
mospheric neutrinos, and in the neutrinos from the sun as well. Since two particle 
cross-sections involving a neutrino and any constituent of the matter are far too small 
to make a substantial effect on the neutrino fluxes, any such effect must come from 
beyond the standard model and also involve some kind of coherence, or other long 
range phenomenon. We consider here the case in which a light scalar particle coupled 
very weakly to both neutrinos and to a constituent of the matter gives rise to neutrino 
masses and mixing within the medium. It is the small mass of the scalar particle that 



1 



allows the very weak coupling to produce the requisite mass differences, while not 
changing two particle cross-sections appreciably. 

We take a real scalar field, 5(x),with mass, /i, interacting with electrons and 
neutrinos through the interaction Lagrangian density, 

Li = GetpelpeS + ^ Gij '. VjViS :, (f) 

where the indices, i,j , run over the three flavors of massless Dirac neutrinos. We 
take the coupling matrix G to be of rank one or two, so that there will always remain 
at least one neutrino combination that is uncoupled to S. We have not prescribed the 
couplings of the scalar field S to fj, , t , since they will not enter into the consider- 
ations that follow. For the case of an otherwise free S field, the interaction (|I|) leads 
in lowest order to a mass matrix, rriij, for neutrinos in a medium of nonrelativistic 
electrons with electron density n e , 

m itj = G e Gijn e fi~ 2 . (2) 

Since the symmetry, v — > 75Z/ , S — > —S, is broken by the coupling to the electrons, 
loops involving electrons will generate a vacuum neutrino mass, beginning in third 
order in the parameter G e . However we shall take G e to be extremely small, and the 
vacuum masses should be negligible. We can choose the couplings and S particle mass 
to make a variety of oscillation-mixing models. We should mention that similar ideas 
have been presented in refs. || and [[jj . In || the scalar field is taken to couple only to 
i/'s, which have a vacuum mass of some other origin, with the possibility of clustering 
to produce effects somewhat like ours. In JF], in which the scalar particle is a dilaton, 
the v mass is given by an expression analogous to (0). Before giving examples in 
our model we consider constraints on the parameters that must be imposed to avoid 
contradictions with what is already known. 

2. S particle luminosity of the sun. Since the S particles will have mass 
of much less than one KeV, they can be produced thermally in the sun; and since 
their interactions with ordinary matter will be far too weak to generate S-opacity 
in the medium, the interaction must be weak enough that the S luminosity of the 
sun, L^ n is a small fraction of the solar constant. We estimate this luminosity by 
considering the free-free, free-bound, and bound-bound processes of S emission and 
absorption, as in the photon case. For both photons and S particles we use the 
lowest term in the expansion in powers of k of the angle averaged (and polarization 
summed, for the photon) squared matrix element between nonrelativistic electronic 
states, where (u, k) is the energy-momentum of the absorbed particle. For a photon 
it is (2/3)e 2 w 2 | < f\r\i > | 2 ; for the scalar particle it is (l/3)G 2 fc 2 | < f\r\i > | 2 . 
Thus for energies of a few KeV, where fcsiw for the S, the inverse absorption lengths 
/I* 1 , and l^ 1 are related by the coupling constant ratio. Detailed balance turns Z" 1 
into a rate of photon energy emission per unit volume. We then scale this rate down 
by the coupling constant ratio. Using median values Z™ 1 ~ .003 cm., T w .8KeV, 
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and radius of 5 x 10 10 cm. [], we estimate a rate of, L^ n ~ 2 x 10 58 G 2 ergs/sec, 
to be compared with the solar constant, 3.90 x 10 33 ergs/sec, leading to the bound, 
G\ << 10~ 25 . This is a smaller bound on G 2 than that found for the axion coupling in 
the estimation of light axion bremsstrahlung luminosities ||, by a factor of about 10 6 , 
for the reason that pseudoscalar coupling in the axion case produces an extra factor 
of roughly T 2 /m 2 e in the rate. Our combined processes have temperature dependence 
roughly proportional to T 1 / 2 ; the weak dependence on temperature makes it unlikely 
that hotter stars will impose more severe restrictions, in contrast to the axion case. 

3. Long range forces For a massless S and a value of G 2 = 10~ 25 the long 
range attraction between lumps of ordinary matter would be roughly 10 12 times that 
of gravity. For fi > 10 _2 eV, or the range of the force less than .002 cm., there appears 
to be no conflict with long range force measurements. However it will turn out that 
these constraints put us on the very border of the domain of values that are needed for 
atmospheric neutrino transformations. f] Defining the model above, with G\ < 10 -25 
and n > 10~ 2 eV, as "weak" , we also consider a less constraining "very weak" form, 
in which G\ < 10 -43 and \i > 10~ n eV; or a range up to 20 km and a strength less 
than 10 -6 that of gravity, putting the parameters well within the bounds allowed by 
the "fifth force" searches in recent years ||. 

4. Atomic physics. The limit on the electron-S coupling in sec. 2 guarantees 
that the forces from S exchange make corrections much less than one part in 10 -23 
to the energy levels of the positronium system and negligible contributions to the 
binding energy per electron of bulk matter. 

5. Supernova neutrino pulse— transparency of the interstellar medium. 
The most stringent bounds on G u are set by the observations of neutrinos from SN 
1987a. It has been shown []n| that the absence of large scattering from relic big bang 
neutrinos (of SN neutrinos coming from the LMC) can be assured with a coupling 
G 2 U < 10 -7 . Thus the upper limit of the combination that enters our v masses is 
G v G e ^r 2 < 10- 10 (eV)- 2 for the "weak" choice, and G u G e fi~ 2 < 10~ 3 (eV)- 2 for the 
"very weak" choice of parameters. Although we choose not to, we could avoid even 
this limitation, since we are taking at least one combination of neutrinos not to couple 
to S. As long as we have a single uncoupled v combination that contains a reasonably 
large fraction of v e the couplings are not limited by the data from 1987a. 

6. Neutrino cross-sections. We demand that the electron-neutrino cross- 
sections produced by S exchange between e_ and v , be significantly smaller than 
standard model cross-sections, for energies down to 1 MeV. Taking G v to be the 



1 We used the opacities plotted as a function of T and p in D. D. Clayton, Principles of Stellar 
Evolution and Nucleosynthesis McGraw-Hill, New York (1968). A multizone estimate using data 
from these figures gives nearly the same limits on coupling. There is a similar scaling of the rate for 
7 + e~ — > S + e~ to the Compton scattering rate in the star, which enables us adequately to take 
into account this process. 

2 The constraints could be relaxed through the agency of an S self- interaction potential, which 
can change the connection between the S field inside a medium and the long range force between 
two separated lumps of matter. 
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maximum eigenvalue of the coupling matrix, and taking only the partial cross-section 
for angles greater than some small minimum gives, 



G e G v < G 2 F {lMeVf « l(T 2i , (3) 

which is satisfied by many orders of magnitude in either of our parameter domains. 
The small angle part gives a correction which is of order m v j \x times the above, which 
again cannot be appreciable. 

7. Supernova neutrino pulse— generation. If we begin with a coupling 
strength and range that gives a mass of 10~ 2 eV for one of the neutrino states in 
ordinary matter, and if in going to superdense matter at, say, 10 12 g/ cm 3 , the mass 
scales up according to ( §), then the neutrino mass in the dense medium would be 
at the GeV level. This cannot be the correct complete description in a domain in 
which the coupling has become effectively so strong, but it does appear that in the 
superdense case one or two of the neutrino combinations would not be available to 
participate in the job of energy and lepton number transport to the surface. However, 
it was shown in [JTTJ] that in the denser parts of the supernova core quite excessive 
neutrino opacities have been used in the calculations of the pulse, due to the neglect 
of interaction effects in the medium. In consequence, the rate of core energy flow of 
the current calculations can be realized with fewer than three neutrino species, with- 
out changing the temperature profile by much. Add to this the ability of the system 
to adapt by changing the temperature profile, and we get a picture with a slightly 
larger and hotter neutrinosphere than in present models []. The temperature would 
decrease more slowly than in the three neutrino model, since there would be only one 
or two neutrino combinations excitable, even at the neutrinosphere densities. As for 
the effect of v — v scattering on the transport of energy out of a neutron star, we 
provisionally agree with the considerations in W\ indicating that it will not make a 



drastic difference. In the region in which the neutrinos move out diffusively we can 
use a frame that is moving at the drift velocity to see that in this frame v-v scattering 
does not change the distribution functions. However, for the larger ranges of our v-S 
coupling the v-v scattering in the interior will create a quasi-equilibrium of otherwise 
sterile UrS, which are effectively required by the v-v interaction to drift at the same 
speed as the v^s. In this case, the pulse would be only 50% non-sterile, at best. If this 
is a problem, it can be completely removed by choosing the Majorana variant which 
we mention later. Another potential problem is an effect on the neutrino spectrum 
as we go from the diffusive region to the free-streaming region. We expect this to 
be small, since an element of neutrino gas does not expand very much in the rather 
narrow transition region. 

8. Coupling schemes. We treat the simulation of the mixing scheme in which 



3 It is larger because at the lower densities of the neutrinosphere the higher energies of the neutri- 
nos emerging from the dense core count more than the interaction-induced reduction in cross-sections 
in determining the interaction length. 
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the atmospheric neutrino data are interpreted through pure — v T oscillations Qand 
the solar neutron problem is addressed through mixing of v e and some other flavor 
(or flavor combination). We take the eigenstates of the matrix G^j to be the following 
combinations 

v\ = cos(a)u e + 2~ 1//2 sin(a)[i/ At + v T ], 
v-i = — sin(a)z/ e + 2 -1 ^ 2 cos(a)[u fl + u T ], 

u, = 2-^[u,-v T ], (4) 

with the respective coupling eigenvalues Gi,2,3 and mass eigenvalues (in the medium) 
m i,2,3 = GeGi,2,3H~ 2 n e . We take G 3 >> G 2 and G\ = 0. Beginning with the 
atmospheric case, we note that the only mixing induced by G3 alone is the mixing 
between \x v and \x T with a maximal mixture, sin(2#) = 1, and a mass difference given 
by 8m pa m 3 , which we shall take in the region, 10 _1 -10 _2 eV. Using mantle electron 
densities in the earth of 10 24 cm -3 in (|2|) this translates into 

G 3 G e {^/leV)- 2 pa[Qx 1(T 12 - 6 x 1(T 13 ], (5) 

which is satisfiable, with rather little leeway for the "weak" model, but with much to 
spare in the "very weak" model. 

The neutrino combination, 1/3, does not enter the solar problem; only the symmet- 
ric combination, — v T mixes with the electron neutrinos. For a particular E v we 
can choose G2 and a to produce a resonance at some particular point in the sun. But 
higher energy z/s will be resonant at higher densities and lower energy z/s at lower 
densities, because of the density dependence of the mass matrix. In a small angle 
MSW model, fits to data are achieved in which the survival probability is almost nil at 
E v = lMeV but increases rapidly as the energy is lowered to the detector threshold, 
because of the resonance region moving into the energy producing region in the sun. 
The survival increases, but more slowly, as the energy is increased from 1 MeV, this 
time because of the loss of adiabaticity. The effect of this loss can be estimated quite 
accurately using the Landau- Zener type calculation of the transition probability be- 



tween the two eigenstates fll5| . In a small angle model in which G2 has been set to 
agree with 5m 2 from MSW at some intermediate density, this transition probability 
for our model is the same function of the resonant density and mass-matrix mixing 
angle as in MSW. But in contrast to the MSW case, it becomes significant at low 
energies, while the restoration of survival at higher energies can be achieved by the 
resonance moving into the energy producing region. It appears that the existing data 



4 Ref. ]13[ discusses the amount of mixing of with v e allowed in three neutrino oscillation 
schemes that fit the SK atmospheric data and are compatible with the CHOOZ data [Q with 
the conclusion that an 18 % admixture in amplitude is allowed. Although we do not explore such 
possibilities, our class of models should allow a considerably larger admixture, since the electron 
density along the path of CHOOZ neutrinos is roughly 1/2 that of the earth's mantle. 
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is less well fit by such a model, with its gentler rise to the left of the minimum and 
sharper rise to the right, than it is by small angle MSW. 

However, in our model a large angle example can provide the following: a survival 
fraction, f s , of about .4, at E=.4 MeV, which increases almost linearly to a value 
f s ~ .5 — .6 at E=15 MeV, depending on choice of parameters, in the absence of 
regeneration in the earth. This is achieved with .62 < sin a < .68 and 4 x 1CT 14 < 
(fj,/leV)~ 2 G2G e < 8 x 10~ 14 . For these parameters, nighttime regeneration can give 
an increase of f s in the region E <2MeV, reaching a 25%-30% increase in f s at .4 
MeV. The all-over behavior is qualitatively consistent with the survival curves for the 
preferred MSW models given in |n|, fig. 9. 

9. Early universe, just before nucleosynthesis. Here the medium consists 
almost entirely of electrons, positrons, neutrinos, antineutrinos, and photons. The 
interaction (|l|) provides long range interactions among all of these species except the 
photons. The induced mass, m^, of the most strongly coupled v combination, coming 
from interaction with the thermal electrons and positrons is, in lowest order, 

m 3 = 4(27r)- 3 G 3 G e m eyU - 2 / = [1 + exp(^ 2 + m^/T)]" 1 , (6) 

Taking G^G e fi~ 2 from (|5|)we find, for example, that in the electron-positron plasma 
at T=2 MeV the range of 7713 is .64-6.4 MeV, so that there would appear to be the 
possibility of largely excluding or largely allowing a sea of u 3 , u 3 during this critical 
period in the expansion. In the early universe there is not sufficient excess of electrons 
over positrons to mix the states (|J) by a significant amount and the mass irt2 induced 
when G3 — > G2 in (j^j will be much smaller. A perturbative estimate of the mass 
changes induced by v — v interactions, through (P) with m e — > m u , will give small 
changes, owing to the neutrino mass factor on the right hand side, f] We conclude 
that there is reduction in the number of effective left-handed v species, most probably 
from 3 to 2, but possibily greater or less. 

There are also new degrees of freedom excited in the soup. S particles equilibrate 
easily from the reaction v + v — * 2S, and right handed z/'s are produced easily in 
v L + u L — > u R + u R , giving effectively, 4.5 species. This can be reduced to a maximum 
of 2.5 by using Majorana fields for the neutrinos, with the scalar coupling the analogue 
of a Majorana mass term. 

10. Conclusion. Until such time as the data rule out our initial qualitative 
assumption, that the downgoing atmospheric i/'s are not altered from the point of 
production, it is worth considering the possibility that neutrinos remain nearly mass- 
less and unmixed in the vacuum. As to the class of models proposed here, the largest 
negative aspect is the introduction of the new particle. The possibility of reducing the 
effective number of neutrinos in the early universe may be a positive aspect, and the 

5 But nonperturbative solutions of the coupled v — v mass equations might lead to much larger 
masses than those discussed above for both v-i and v%. 



6 



extension of the duration of the supernova pulse, without diminishing the effective 
temperature, may be one as well. 
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